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Mixing and transport of passive particles are studied in a simple kinematic model of a mean- 
dering jet flow motivated by the problem of lateral mixing and transport in the Gulf Stream. We 
briefly discuss a model streamfunction, Hamiltonian advection equations, stationary points, and 
bifurcations. The phase portrait of the chosen model flow in the moving reference frame consists of 
a central eastward jet, chains of northern and southern circulations, and peripheral westward cur- 
rents. Under a periodic perturbation of the meander's amplitude, the topology of the phase space 
is complicated by the presence of chaotic layers and chains of oscillatory and ballistic islands with 
sticky boundaries immersed into a stochastic sea. Typical chaotic trajectories of advected particles 
are shown to demonstrate a complicated behavior with long flights in both the directions of motion 
intermittent with trapping in the circulation cells being stuck to the boundaries of vortex cores and 
resonant islands. Transport is asymmetric in the sense that mixing between the circulations and 
the peripheral currents is, in general, different from mixing between the circulations and the jet. 
The transport properties are characterized by probability distribution functions (PDFs) of durations 
and lengths of flights. Both the PDFs exhibit at their tails power-law decay with different values of 
exponents. 

PACS numbers: 05.45.Ac; 47.52. +j; 92.10.Ty 



I. INTRODUCTION 



Passive particle advection can be considered as a Hamiltonian dynamical problem. Advection in 
a two-dimensional meandering shear flow, a simple kinematic model of lateral mixing and transport 
in the Gulf Stream, is an example of a system with 3/2 degrees of freedom. The phase space of an 
advected particle in such a flow is mixed, i.e. it consists of chaotic layers, islands of periodic motion, 
and layers of periodic open trajectories. All these parts of the phase space influence the transport of 
passive particles, which can become anomalous if trajectories stick to the boundaries of topological 
structures with regular dynamics. Particularly the stickiness to the so-called ballistic islands can create 
flights along the stream flow, make the transport superdiffusive, and influence in a strong way the 
advected particles distribution. 
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The goal of this paper is to study Lagrangian lateral mixing and transport of passive particles advected by a shear 
flow that has a coherent space organization with two chains of vortices and a strong jet between them. Temporal 
perturbations of shear flows are known to induce chaotic mixing. Similar flows occur naturally in the ocean as western 
boundary currents, like the Gulf Stream in the Atlantic Ocean and the Kuroshio in the Pacific Ocean, separating 
water masses with different physical and biogeochemical characteristics. 

Hamiltonian chaos theory has been widely used to study advection of passive particles in hydrodynamic flows 
of ideal fluid (for review see 

HQ). 

The equations of motion of a passive particle advected by a two dimensional 
incompressible flow is known to have a Hamiltonian form 

dx . .9* 
— = u(x,y,t) = - — , 

dt dy m 

dy 9M> K 1 

- =v{x , y>t ) = —, 

with the stream function \& playing the role of a Hamiltonian and the coordinates (x, y) of a particle playing the 
role of canonically conjugated variables. Equations form a Hamiltonian dynamical system whose phase space is 
identified as the physical space of advected particles. 

All time-independent one-degree-of-freedom Hamiltonian systems are known to be integrable. It means that all 
fluid particles move along streamlines of a time-independent streamfunction in a regular way. Equations with a 
time-periodic streamfunction are usually non-integrable, giving rise to chaotic particle's trajectories. In other words, 
even regular Eulerian velocity fields (u, v) can generate chaotic trajectories, the phenomenon known under the names 
"chaotic advection" or "Lagrangian chaos" In this field mixing means stretching and deformation of material 

lines by advection. The term "transport" refers to the motion of passive particles from one region of the physical 
space to another. 

Well developed theory of Hamiltonian chaotic mixing and transport in the phase space (for a review see 
has been successfully applied to study chaotic mixing and transport of passive particles in fluids lull 
ll-U. 11.4 lid Il7t ll8L llfll l2("l | . The phase space of a typical chaotic Hamiltonian system consists of regions of chaotic 
mixing and regions with regular trajectories. A boundary between them often contains a complicated hierarchical 
structure of islands with regular trajectories and can be extremely sticky for some values of control parameters 
[3, 13 0| • All these properties are manifested in mixing and transport of passive particles in fluids and 
have been observed in laboratory [23, 12a, l22( • Particles in the core regions of regular motion and in the islands are 
trapped there forever unless we do not take into account molecular diffusion and noise. They do not participate in 
the transport. So the outermost KAM tori play the role of transport barriers in the fluid. The particles, entering the 
boundary layer, can spend a very large time there giving rise anomalous transport or anomalous diffusion. 

The motion in chaotic regions is extremely sensitive to small variations in initial particle's positions. So one forces 
to use an statistical approach to describe transport. A commonly used statistical measure of transport is the variance 
0-2 (t) = ( x2 ), which for simplicity is written for the particle displacement in one dimension. The averaging is supposed 
to be done over an ensemble of particles. If mixing in the phase space would be ideal (as it occurs in a hyperbolic 
dynamical system with unstable orbits only, whose phase space does not contain any regular islands) the variance 
would satisfy the law of normal diffusion, u 2 {t) = 2Dt, with a well-defined diffusion coefficient D — lim a 2 /2t. If 

t — >oc 

trajectories are dominated by sticking to the boundaries of regular regions, where particles may spend a long time, 
subdiffusion takes place with the transport law a (t) ~ f, /i < 1. Superdiffusion with fi > 1 occurs when particles 
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execute long "flights", i.e. travel long distances between (or even without) sticking events. Anomalous diffusion 
implies the diffusion coefficient D to be either zero (// < 1) or infinite (fx > 1). As to transport of advected particles 
in fluids with diffusion Q and transport in the phase space in the presence of noise [sil l^. the situation is more 
compicated. For example, as it is shown in Ref. [^ij subdiffusiion is not possible in incompressible flows in the 
presence of molecular diffusitivity. 

It is as well possible to describe transport in terms of spatial extensions and durations of sticking and/or flight 
events. The respective probability density functions (PDFs) of durations of flights, P(t) ~ £~ 7 , and of lengths 
of flights, P{x) ~ x~ v , are expected to be power-law functions. Long-time flights provide a link between chaotic 
advection and anomalous diffusion. Quantitative connections between all the exponents /i, 7, and v have been 
established with some simple models of chaotic motion and continuous-time random walks (for a review of anomalous 
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transport see |9.l33U34l35l|). 

To specify the model we choose the shear flow in the form of the Bickley jet {jjfj with the velocity profile sech 2 y 
and an imposed running wave with the periodically varied amplitude. This model is related to the so-called non-twist 
Hamiltonians (see [L"! I37I and references therein) that have a specific phase space topology and bifurcations with 
respect to the parameters change. As it usually occurs in Hamiltonian systems, separatrices are destroyed in the 
presence of an arbitrary small time-dependent perturbation. Instead of them narrow stochastic layers appear which 
are "seeds" of Hamiltonian chaos. There are an infinite number of nonlinear resonances between the frequency of 
the perturbation and frequencies of rotations of advected particles. With increasing the strength of the perturbation, 
more and more of them begin to destroy and overlap. Motion of particles in the regions of the overlapping resonances 
in the phase space is no more confined by impenetrable barriers of KAM tori and a large-scale transport becomes 
possible. 

In Sec. m we derive equations of motion for passive particles in the moving frame of reference, find stationary points 
of the flow, topologically different regimes of motion, and bifurcations between them. Similar kinematic models with 
imposed time-periodic velocity fields have been used in literature on physical oceanography to study lateral Lagrangian 
mixing and transport in the meandering Gulf Stream current 0, Q, Q] ■ Mixing has been shown with the help of 
Poincare sections to be chaotic in a wide range of the control parameters, and the magnitudes of fluxes have been 
estimated with the help of the Melnikov integral. 

In Sec. [TO] we describe topology of the phase space with regions of regular and chaotic motion of advected particles. 
Sticking to the boundaries of different types of islands, including ballistic ones, and of the vortex cores is demonstrated 
on Poincare sections of the particle's trajectories. Sec. |^ is devoted to quantifying chaotic transport. We compute 
the statistics of durations and lengths of flights for a number of long-time trajectories and find the transport both in 
western and eastern directions to be anomalous with different values of the exponents. 
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II. MODEL STREAMFUNCTION, STATIONARY POINTS, AND BIFURCATIONS 

We take the Bickley jet with a running wave imposed as a model of a meandering shear flow in the ocean Q. The 
respective streamfunction in the laboratory frame of reference has the following form: 



y'{x',y',t') = -*^tanh 



y' — a cos k(x' — ct') 
v xJl + k 2 a 2 sin 2 k{x' - ct') ( 



(2) 



where the hyperbolic tangent produces the velocity profile ~ sech 2 y' , the square root provides a possibility to have 
a homogeneous Bickley jet with the width A, and a, k, and c are amplitude, wavenumber, and phase velocity of the 
wave, respectively. The normalized streamfunction in the frame moving with the phase velocity c is 



, / y — A cos x \ „ 
* = - tanh — -. + Cy, 

\L^l + A 2 sm 2 xj 
where x = k(x' — ct'), y — ky' , A — ak, L — Xk, and C = c/ty' k. 



(3) 




FIG. 1: The phase portrait of the model shear flow 10 in the frame moving with the meander's phase velocity. Streamlines 
in the circulation (C), jet (J), and peripheral currents (P) zones are shown. The parameters are A = 0.785, C = 0.1168, and 
L = 0.628. 

The respective advection equations Q in this moving frame have the form 

-c, 



+ A 2 sin 2 a; cosh 2 
A sin x(l + A 2 — Ay cos a;) 

V ~ ~L(l + A 2 sm 2 x) 3/2 cosh 2 6 (4) 

y — A cos x 

9 = — 7 

LV 1 + A 2 sin 2 x 

where dot denotes differentiation with respect to dimensionless time t = ^' Q k 2 t' . The one-degree-of-freedom dynamical 
system (2J is generated by the Hamiltonian and has three control parameters — the jet's width L, meander's 
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amplitude A, and phase velocity C. The scaling chosen results in the translational invariance of equations (@J along 
the x-axis with the period 2n. 

A brief description of stability of stationary points and bifurcations of equations <Q} is given below. It is clear from 
the first equation in iQJ that stationary points may exist only if the condition LC < 1 is fulfilled. Two equalities, 
following from the second equation in <QJ 



sinx = 0, 1 + A 2 - Ay cos x = 0, (5) 



give four stationary points 



#1,2 = 0, yi 2 = ±L Arcosh J + A, 

V i£ (6) 

x 3A = 7r, 2/3,4 = ±L Arcosh J — - A. 

V i^O 

As follows from the stability analysis, the second and third points are always stable whereas the first and fourth ones 
are stable only if 

AL Arcosh \ — > 1. (7) 
V IjLj 



C < 5 (8) 

Lcosh 2 (l/AL) v ; 



It is proved in Ref. |42j that if the condition 



is fulfilled, then there are four new stationary points in equations iQJ in addition to the points 10. 

Having found the stationary points, their stability properties, and how they are connected, we get the following 
topologically different phase portraits of the streamfunction Q 

1. C > C cr i = l/L, there are no stationary points. 

2a. C cr i > C > C CI 2 = l/£cosh 2 (l/AL) and C > C cr 3, there are four stationary points two centers (the second 
and third ones) and two saddles (the first and fourth ones). There are two separatrices each of which passes 
through its own saddle. A free flow between the separatrices is westward. 

2b. C cr i > C > C cr 2 and C < C cr 3, the stationary points are the same as in the preceding case, but a free flow 
between the respective separatrices is eastward. 

3a. C cr 2 > C > C cr 3, there are eight stationary points: four centers © and four saddles with positions found 
in Ref. ^^|. There are two separatrices, each of which connects two saddle points. A free flow between the 
separatrices is westward. 

3b. C cr 2 > C and C < C cr 3, the stationary points are the same as in the preceding case, but a free flow between 
the respective separatrices is eastward. 

Therefore, from the point of view of existence and stability of stationary states, there are three possibilities: 1) 
there are no stationary points; 2) there are four stationary points, two centers and two saddles; 3) there are eight 
stationary points, four centers and four saddles. A bifurcation between the first and second regimes consists in arising 
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two pairs "saddle-center". A bifurcation between the second and third regimes consists in arising two saddles and a 
center between them instead of one saddle (a fork-type bifurcation) . 

There is one more bifurcation that does not change the number and stability of stationary points but changes the 
topology of the flow. The values of the streamfunction (01 on the separatrices are equal on modulo but of opposite 
signs. There is a critical value of the phase velocity C = C cr3 under which the separatrices coincide and the respective 
streamfunction is equal to zero. If C > C = C cr 3, a free flow between the separatrices is westward, whereas with 
C < C = C cr3 it is eastward. It is difficult to find C = C cr3 analytically but it may be shown 0] that C cr3 > C cr2 , if 

2(1 + A 2 ^ 



Otherwise C cr3 < C, 



AL sinh(2/AL) 



< 1. 



(9) 



^cr2- 




FIG. 2: The frequency map represents by color the values of frequency of rotation / of particles with initial positions (xo, Do) 
advected by the unperturbed flow Q. 



III. TOPOLOGY OF THE PHASE SPACE 

Being motivated by the boundary oceanic eastward meandering jet currents like the Gulf Stream and the Kuroshio, 
we will deal further in this paper with the case 2b in the list in the preceding section. The respective phase portrait 
in the reference frame, moving with the phase velocity of the meander (Fig.EJ, consists of three different regions: the 
central eastward jet (J), chains of the northern and southern circulations (C), and the peripheral westward currents 
(P). The elliptic points (centers) of the circulations are at critical lines y — y c with u(y c ) = c and v(y c ) — 0. The 
northern separatrix connects the saddle points at x\ = 2irk and y\ to be defined from equation @ and the southern 
one connects the saddle points at xa = [2k + l)ir and j/4 from (O, where k = 0, ±1, .... 

In the presence of a time dependent perturbation, the separatrices are destroyed and transformed into stochastic 
layers. As a perturbation, we take the periodic modulation of the meander's amplitude, A = Aq + e cos(ujt + ip) , in the 
streamfunction © written in the frame moving with the phase velocity of the meander. The jet's width L — 0.628, 
the phase velocity C = 0.1168, and the meander's amplitude Aq = 0.785 will be fixed in simulation throughout the 
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FIG. 3: (color online). Poincare section of the streamfunction 10 with the same parameters as in Fig. U and the amplitude 
£ = 0.0785 and frequency of perturbation w — 0.2536. (a) General view, (b) zoom of a pair of islands of the secondary resonance. 

paper. These values are in the range of the parameters which has been estimated 0,13 to be realistic with the Gulf 
Stream current. 

Changing the strength e and frequency to of the perturbation, we can change the number of overlapping resonances 
between the perturbation frequency and the frequencies / of particle's rotations in the circulation zones. In Fig.[2]we 
plot a frequency map f(xo,ya) that shows by color the values of / for particles with initial positions (xo,yo) in the 
unperturbed system. 

Let us take the values of the perturbation frequency to = 0.2536 to be close to the values of the unperturbed frequency 
/ in the inner core of the circulation zone and the strength of perturbation to be rather small e — 0.0785. Due to 
the zonal and meridian symmetries it is possible to consider particle's motion on the cylinder with < x < 2tt. The 
respective Poincare section for a large number of trajectories is shown in Fig. 03 The vortex core, that survives under 
this small perturbation, is immersed into a stochastic sea where one can see six islands of a secondary resonance to be 
emerged from three islands of the primary resonance 3/ = 2io, where / = 0.169. A pair of the secondary resonance 
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FIG. 4: A single ballistic trajectory (a) sticking to the vortex core (b). Bold dots show particle's positions through the period 
of the perturbation. 



islands is zoomed in Fig. Et>- 

Sticking of particles to boundaries between regular and chaotic regions in the phase space is a typical phenomena 
with chaotic Hamiltonian systems Q, B, 13 13 13 Q] ■ What is special in chaotic advection, that real fluid parcels 
may be trapped for a long time near a vortex core and borders of resonant islands. As to the vortex core, we illustrate 
this phenomenon in Fig. 01 where tracks of a single ballistic trajectory (Fig.QJi) through the period of the perturbation 
2tt/uj are shown by bold dots (Fig. \%}p). In Fig. we plot tracks of another ballistic trajectory (Fig. Et) sticking to 
the islands of the secondary resonance. 

Besides the resonant islands with particles moving around the elliptic point in the same frame (Fig. , we have 
found so-called ballistic islands situated in the chaotic sea, in the peripheral currents and on the border between the 



-0.5 



12 



18 



24 



1.06 



1.04 



1.02 



..•'1 *.j 
" ■ '< -.V 

■ " " ■*•■ 



1.97 



1.98 



1.99 



FIG. 5: A single ballistic trajectory (a) sticking to two of the resonant islands of the secondary resonance (b) Bold dots show 
particle's positions through the period of the perturbation. 

chaotic sea and the meandering jet. Ballistic modes correspond to the stable periodic motion of particles 

from one frame to another. When mapping their positions at the moments of time tk — 2kir /u) (k=l,2,. . . ) onto the 
first frame, one can see chains of ballistic islands along the stochastic layer that replaces the destroyed separatrix 
(Fig.Et). Zoom of one of the ballistic islands and stickiness of a chaotic trajectory to its border are shown in Fig.^b- 
In Fig. [3 we demonstrate two different ballistic trajectories. If a particle is placed initially inside the ballistic island, 
it travels (to the west in this case) in a regular way (see the upper trajectory in Fig. Ek) with practically a constant 
or slightly modulated zonal velocity (see the lower line in Fig.0}). The lower trajectory in Fig. Ct- on the (x, y) plane 
corresponds to a particle placed initially nearby the border of the same ballistic island from the outside. Dependence 
of its zonal position on time (the upper curve in Fig. [3d) demonstrates clearly an intermittency of flight and sticking 
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FIG. 6: (color online), (a) Chains of ballistic islands on both sides of the border between the peripheral current and the 
circulation zone, (b) Zoom of one of the ballistic islands and its sticky border. 

events. Stickiness to the islands in Figs. El - El is shown for different trajectories, for a convenience. Actually, any 
chaotic trajectory sticks to different islands at different time intervals. Such a behavior one can see in Fig. This 
type of dynamics complicates its description because of a multifractal nature. 

IV. CHAOTIC TRANSPORT 

With increasing the perturbation strength e, the boundaries (see Fig.QJ between the jet (J), circulations (C), and 
peripheral currents (P) begin to break down. The vortex core shrinks (see Fig. EJ, oscillatory and ballistic islands 
appear (see Figs. and EJ, and a stochastic layer arises instead of the impenetrable transport barrier, the unperturbed 
separatrix. Passive particles which stick to the boundaries of regular motion can spend there a significant time and 
deposit to subdiffusion. Those which find themselves in the jet may travel long distances to the east between sticking 
events, whereas the particles in the peripheral currents travel in the west directions. Because of the absence of the 
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PIG. 7: (color online). Examples of ballistic trajectories (a) on the {x,y) plane and (b) on the (t,x) plane. The blue regular 
trajectory, which is the upper one in (a) and the lower one in (b), is inside a ballistic island. The green chaotic trajectory 
with close initial position, which is the lower one in (a) and the upper one in (b), demonstrates intermittent sticking and flight 
events. 



former transport barriers, the motion of some particles is intermittent with a large number of turning points at which 
particles change directions. As usual, we will characterize statistical properties of particles motion by probability 
distribution functions (PDFs). We will call "a flight" any event between two successive changes of signs of the 
particle's zonal velocity. In this terminology a sticking consists of a number of flights with approximately equal flight 
times. It should be noted that some authors distinguish flights from sticking motion by examining a distance between 
successive local extrema of trajectories x(t) QjQl- Both the PDF of durations of flights P(Tf) and of lengths of 
flights P{xf) will be used to characterize the chaotic transport. 

The Poincare section of the flow with the perturbation strength e = 0.0785 and frequency ui = 0.2536 is shown in 
Fig. Et. The flight PDFs are computed with 6 particles (initially placed in the first east frame inside a stochastic 
layer) up to the long time t = 5 ■ 10 8 and with 125 particles up to t — 5 • 10 6 . Particles inside the stochastic layer 
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FIG. 8: (a) Probability distribution functions (PDFs) of lengths of eastward (triangles) and westward (squares) flights obtained 
with 6 chaotic trajectories for 5- 10 s steps, (b) For short eastward flights, 15 < Xf < 85, the PDF decays exponentially, whereas 
the PDF for westward flights oscillates in the same range, (c) The PDF tail of lengths of westward (squares) long flights decays 
as a power law with the slope v = 3.12 ± 0.1. 

may execute very complicated motion changing the direction in a chaotic way, sticking to the island's boundaries and 
executing ballistic flights with different values of length and duration. 

The PDF of the lengths of flights is shown in Fig. [Hk for both the directions. The asymmetry between the eastward 
(shown by triangles) and westward (shown by squares) flights is evident. The transport velocity in the peripheral 
currents is different from the transport velocity in the jet. Mixing between the circulations and the peripheral currents 
is, in general, different from mixing between the circulations and the jet. Both the PDFs can be roughly split into 
three distinctive regions. The very short flights with small values of \xf\ (< 2ir) are supposed to be dominated by 
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FIG. 9: The number of particles N executing short flights with \xf\ < 2tv and a fixed duration Tf. Peaks at Tf ~ 11.8 and 
Tf ~ 18 correspond to particles sticking to the boundaries of the vortex core and oscillatory islands. 

sticking to the boundaries of the vortex core and oscillatory islands. To check that we plot in Fig. El the number of 
particles N executing short flights with small values of Tf and with the lengths of flights \xf\ < 2ir. The prominent 
peaks at Tf ~ 11.8 and Tf ~ 18 correspond to those particles which stick to the vortex core and resonant island's 
boundaries (see Fig. Eland QJ. The PDF for eastward flights with the lengths in the range 15 < Xf < 85 decays 
exponentially, whereas the PDF for westward flights is an oscillating function in this range (see Fig.|Ht>). The tails 
of both the PDFs are close to a power-law decay Pf(x) ~ | | ^ . In Fig. |Ht we estimate the exponent for long 
westward flights with corresponding error by least-square fitting of the straight line to the log-log plot of the data 
to be v = 3.12 ± 0.1. The tail for long eastward flights in the same range is an oscillating function whose slope is 
difficult to measure because of insufficient statistics. Simulation with 125 initial conditions up to the time t — 5 • 10 6 
shows practically the same results. 

The PDF of durations Tf of eastward flights is shown in Fig. with different number of trajectories to be 
computed and up to different times. In Fig.llOfa, we plot the time PDF with 125 chaotic trajectories for 5 • 10 6 steps. 
The exponent of the function Pf(t) ~ i~ 7 is estimated to be 7 = 2.43 ± 0.1 in the range from log 10 Tf = 2.75 to 
log 10 Ty = 3.8. In Fig.llOb the time PDF is plotted with 6 chaotic trajectories but for much longer times 5 • 10 8 . The 
respective exponent is estimated to be practically the same, 7 = 2.48 ± 0.05. What is remarkable is that both the 
time PDFs demonstrate at the very tails periodic oscillations in the logarithmic scale. That is a common feature with 
time distribution functions in typical Hamiltonian systems . One of the possible explanations of this fact is 

existence of a fractal-like hierarchy of islands [6| and a discrete renormalization group j3| ■ The respective PDF of 
durations of westward flights is shown in Fig. El with 125 chaotic trajectories for 5 ■ 10 6 steps (Fig. \llh ) and with 
6 chaotic trajectories for 5 • 10 s steps fFig. ITThV The single algebraic tail is not so evident with westward flights as 
with eastern ones. Rather one can see in this case more prominent log-periodic oscillations. 
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FIG. 10: The PDF of durations of eastward flights (a) with 125 chaotic trajectories for 5 • 10 6 steps (the slope is 7 = 2.43 ±0.1) 
and (b) with 6 chaotic trajectories for 5 • 10 8 steps (the slope is 7 = 2.48 ± 0.05). 

As it is seen from Fig.^3 the flight-time distribution Pf(t) has the asymptotic 

, , const , 
Pf(t) = (10) 

One can consider it as an intermediate asymptotic depending on which set of sticky islands it represents. We can 
expect a fairly complicate multifractal dependence of Pf on t and the simplified formula itTflll should be applied to 
specific interval of time and parameters (see more discussion in Refs. 

HQ) 

. Another approach to study the stickiness 
is to consider a distribution of Poincare recurrences of trajectories P lec (t, B) to a small domain B taken in a stochastic 
sea far enough from sticky domains. It was shown in Refs. (fl,LD| that for many different models of chaotic Hamiltonian 
dynamics the stickiness leads to a distribution of density probability to find the recurrence time within the interval 

(t,t + dt): 

const 

P™(t,B) = — sr , (11) 

where 7 rcc = min7, (t —* 00). From Fig. we have 7 = 2.43 ± 0.1) and 7 = 2.48 ± 0.05). This result is consistent 
with the Kac lemma that imposes 7 > 2. 

A typical chaotic particle in a stochastic layer changes many times its direction (the sign of the zonal velocity u) . 
We have found that the respective turning points are in a narrow strip confined by the curves to be defined by the 
unperturbed equations of motion J3J with x = and A = Aq ± e. As an example, we demonstrate in Fig. I12h (color 
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FIG. 11: The PDF of durations of westward flights (a) with 125 chaotic trajectories for 5 • 10 6 steps and (b) with 6 chaotic 
trajectories for 5 ■ 10 s steps. 



online) a narrow strip of red turning points confined by the curves with Aq = 0.785 and e — 0.0785. All these points 
belong to a single chaotic trajectory on the cylinder < x < 2n. Figure fl2h provides a zoom of the lower part of the 
strip with the confining curves 



= lVi + a 2 sin 2 x Arcosh , 



sin 2 x 



A cos : 



(12) 



and A — Aa ± e to be shown by black solid lines and three regular trajectories to be shown by different colors. The 
empty part of the strip corresponds to the place where the vortex core is situated (see Fig. Et) • 



V. CONCLUSION 

In this paper we have investigated mixing and transport of passive particles in a two-dimensional nonstationary 
shear flow of the ideal fluid that is relevant to study of mixing of water masses along with their properties in jet-like 
western boundary currents in the oceans. 

We have considered a known kinematic model of a meandering sech 2 zonal flow 0,13' derived the respective 
advection equations in the reference frame, moving with the phase velocity of the meander, found the stationary 
points and their stability, and possible bifurcations. The control parameters of the model streamfunction the 
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FIG. 12: (color online). Turning points (red ones online) of a single chaotic trajectory on the cylinder < x < 2ir are in a 
narrow strip confined by two curves with equations y = f(x) to be defined in the text, (a) General view, (b) Zoom with a few 
regular trajectories shown for reference. 

scaled jet's width L = 0.628, the phase velocity C — 0.1168, and the meander's amplitude A = 0.785, were chosen to 
correspond to estimated realistic values of the Gulf Stream meander's parameters 0,0] and to the topology of the 
flow with two chains of circulation cells separated by an eastward jet (the case 2b in the list of possible flow regimes 
in SecHD. 

The respective phase space has been shown to consist of the vortex cores filled with regular trajectories surrounded 
by a chaotic sea with stochastic layers and chains of islands of oscillatory and ballistic regular motion. The boundaries 
of these regular structures have been shown to be sticky. Typical adverted particles may alternate chaotically between 
being trapped near the vortex cores and island's boundaries and "flying" freely in the jet. 

Defining "a flight" as an event between two successive changes of signs of the particle's zonal velocity, we characterize 
the statistics of flights by the PDFs computing both the PDF of durations of the flights Pf(t) and of lengths of the 
flights Pf{x). Both the functions have been found to have regions with different laws. For short flights we have found 
two peaks in the Pf{x) function corresponding to sticking events in the circulations, for middle flights an exponential 
decay has been found, and the tails with long flights has been found to fit to a power law. An asymmetry of transport 
of passive particles in the eastern and western directions is caused by a difference in mixing between the circulations 
and the peripheral currents and mixing between the circulations and the jet. The tails of the time PDFs for eastern 
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and western flights have been shown to be close to power-law decay functions. The anomalous statistics of the flights 
is caused by the coherent structures - the circulation cells, the peripheral currents, and the jet - where particles can 
spend an anomalous amount of time either being trapped in the circulations or moving in the peripheral currents and 
in the jet. 
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